d-wave to s-wave to normal metal transitions in disordered superconductors 
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We study suppression of superconductivity by disorder in d-wave superconductors, and predict 
the existence of (at least) two sequential low temperature transitions as a function of increasing 
disorder: a rf-wave to s-wave, and then an s-wave to metal transition. This is a universal property 
of the system which is independent of the sign of the interaction constant in the s-channel. 



PACS numbers: 74.20.Mn, 74.72.-h, 79.60.-i 

Generally the order parameter in superconductors is a 
function of two coordinates and two spin indices. Classi- 
fication of possible superconducting phases in crystalline 
materials was done in P, A majority of low- Tc crys- 
talline superconductors have a singlet order parameter 
with s-wave symmetry. It does not change its sign under 
rotation, and in the isotropic case can be approximated 
by a complex number A*(r) — A(r, r). However, over the 
last decades a number of superconductors have been dis- 
covered in which the order parameter changes sign under 
rotation. A notable example is HTC superconductors, 
where in the absence of disorder the order parameter has 
singlet d-wave symmetry 0, Q: A(r - r') chan ges sign 
under rotation by tt/2, and consequently A(r,r) — 0. 
This means that the Fourier transform A(k) changes its 
sign under a n/2 rotation as well, as is shown schemat- 
ically by the rosettes in Fig.l. Since the sign of A(k) 
in crystalline rf-wave superconductors depends on the di- 
rection of the wave vector k, they are much more sen- 
sitive to disorder than s-wave superconductors: at tem- 
perature T = 0, d-wave superconductivity gets destroyed 
when the electron mean free path I is of the order of the 
zero temperature coherence length in a pure supercon- 
ductor, I ^ Iq = 1.78^o 3> 1/kp. Here kp is the Fermi 
wavelength. This is in contrast with the case of s-wave 
superconductors, where according to the Anderson the- 
orem the superconductivity is destroyed at much higher 
level of disorder, when I ~ 1/kp. The fate of the d-wave 
superconductors at ^ < depends on the sign of the in- 
teraction constant As in the s-wave channel. If the inter- 
action As in the s-wave channel is attractive, but weaker 
than the attraction in the d-wave channel |As| < |Ad|, 
then at weak disorder, (I > i^o), the superconducting or- 
der parameter has d-wave symmetry, while at Z < the 
disorder destroys the d-wave superconductivity and the 
system undergoes a phase transition into an s-wave su- 
perconducting state. (See, for example, Q). 

In this article we consider a more interesting case, in 
which the interaction in the s-channcl is repulsive at 
strong enough disorder 1/kp I ^ £,o the system is 
in normal state. We predict at least two low- temper- 



ature phase transitions: a d-wave to s-wave, and then 
an s-wave to normal metal transition. Qualitatively the 
phase diagram of disordered d-wave superconductors is 
shown in Fig.l. Let us first discuss the definition of s- 
and d- symmetries in bulk disordered systems. Before 
averaging over random realizations of disorder, the sys- 
tem does not possess any particular spacial symmetry at 
all. However in bulk samples, the symmetry is restored 
upon configuration averaging. We can think of several 
different definitions of the global symmetry of the or- 
der parameter: a) An operational definition is provided 
by the result of a phase sensitive experiment, such as 
the corner S QUID e xperiment, for example, [3, b) 
The quantity A(r, r') can be characterized as having d- 
wave or s-wave symmetry. Here the over-line stands for 
the averaging over the sample volume, c) A globally s- 
wave component of the order parameter can be defined in 
terms of the local s-component of the anomalous Green 
function jr(r = r') = JF(*)(r). If we define P± to be the 
volume fraction of a sample where F^^\r) has a positive 
or negative sign, respectively, then the system has an 
s-wave component if {P~^ — P~) ^ 0. These definitions 
may be not equivalent under all circumstances. However, 
for the most part, we will deal with the interval of pa- 
rameters in which all these definitions are approximately 
interchangeable. 

It is important to realize that it is inevitable near crit- 
icality to have a situation in which the local pairing in 
disordered superconductors is "d- wavelike" and yet the 
global superconductivity has s-wave symmetry. The d- 
wave to s-wave transition can be understood at the mean 
field level. The electron mean free path is an average 
characteristic of disorder. Let us introduce a "local" 
value of the mean free path l{r) averaged over a size of 
order ^o- In the region of parameters where d-wave su- 
perconductivity is sufhciently suppressed by disorder, the 
spatial dependence of the order parameter can be visual- 
ized as a system of superconducting puddles with anoma- 
lously large values of the order parameter, which are con- 
nected by Joshepson links through non-superconducting 
metal. The superconductivity inside the puddles may 
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be enhanced because either the electron interaction con- 
stant, or the mean free path in the puddles (or both) may 
be larger than their average values. 

Let us assume that the distance between the puddles 
is larger than both their size and the mean free path. In 
this case the system is already in a state with the " global 
s-wave" symmetry. Its origin is illustrated qualitatively 
in Fig. 2, where a system of superconduting puddles of 
arbitrary shape embedded into a metal is shown. The 
order parameter inside the puddles has rf-wave symme- 
try, and the orientation of the gap nodes is assumed to 
be pinned by the crystalline anisotropy. In a rf-wave su- 
perconductor, in addition to an overall phase of the order 
parameter, there is an arbitrary sign associated with the 
internal structure of the pair wave function. Specifically, 
we adopt a uniform phase convention such that when the 
phase of the order parameter (pi — 0, this implies A(r, r') 
in puddle i is real and has its positive lobes along the y 
axis and its negative lobes along the x axis. 

The inter-puddle Joshepson coupling originates from 
the proximity effect in the normal metal. It is char- 
acterized by the anomalous Green function ^(r, r') = 
F{r,r',t = t'), which is connected to A(r,r') by the 
interaction constant. Due to the lack of symmetry at 
the boundary of a puddle, an s-wave component J-{r = 
r') = jr(*)(r) 7^ of the anomalous Green function is 
generated in the neighboring metal. At a distance from 
the superconductor-normal metal boundary larger than 
the elastic electron mean free path the anomalous Green 
function becomes isotropic. In other words, only the s- 
component J-{r — r') = JF^*) (r) survives. It is this com- 
ponent that propagates between far separated puddles 
and determines the Joshepson coupling. 

The sign of .7^*^"^ (r) at a normal metal-superconductor 
boundary, is determined by the sign of the d-wave order 
parameter in the k-direction perpendicular to the bound- 
ary. Therefore it changes along the boundary of a puddle. 

At a distance from an individual i-th puddle larger 
than its size and smaller than the distance between the 
puddles the quantity J-''^{r) has a sign rji = ±1, which 
depends on the shape of the i-th puddle. This point is 
illustrated in Fig. 2a, where the sign of the anomalous 
Green function is positive in hatch-marked areas, and 
negative outside of these areas. 

If the distance between puddles is larger than their 
size, the sign of the Joshepson coupling energy Ejos is 
determined by a product rjirjj, 
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Here indexes i,j label puddles, J^p > 0. Eq. [T] repre- 
sents the Mattis model, which is well known in the theory 
of spin glasses [llj . The ground state of this model cor- 
responds to 



Thus the distribution of cos{(j)i) between puddles looks 
completely random as it is shown in Fig. 2a. However 
the system is not a glass because it's ground state has 
a hidden symmetry. In other words if the distances be- 
tween puddles are bigger than the characteristic size of 
the puddles, R, the Josephson coupling between puddles 
inevitably favors globally s-wave superconductivity, even 
though the order parameter on each puddle looks locally 
rf-wave -like. It is obvious that at a high concentration 
of puddles, the order parameter in the ground state has 
global rf-wave symmetry (See Fig. 2b.). 

At intermediate distances, the situation is more com- 
plicated. Areas with different signs of JF'^''(r) mix in a 
random fashion. We argue that the most important as- 
pects of this complex situation can be modelled by adding 
to the right hand side of Eq. [1] a term 
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where jj^j ^ > characterizes the strength of the exchange 
interaction between the d-wave components of the order 
parameter. Typically, at small jr^ — r^l, J^^^ > J^^\ but 
at large jr^ — r^l the coupling strength J^J-* decays more 

slowly than J^j'' ■ Here are coordinates of the puddles. 
Thus it is likely that in this intermediate region the sys- 
tem may exhibit spin glass features and/or coexistence 
of d-wave and s-wave ordering. In this article, however, 
we will not further explore this fascinating but complex 
aspect of this problem. 

To quantify the picture presented above one has 
to compute the Josephson coupling between a pair 
of far separated puddles. Since the time that it 
takes for electrons to travel between puddles is shorter 
than the characteristic time of fluctuations of the or- 
der parameter on individual puddles, one can cal- 

(s) 

culate using the mean-field Usadel equation for 

the configuration-averaged anomalous Green function 
{Ft ''\r)) = — «sin6'(e,r) in the metal. 



(4) 
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Here Dtr is the transport diffusion coefficient of elec- 
trons in the metal, J-e^\r) is the Fourier transform of 
F^''^[r,{t - t'))], A(^) = A(^)jP'(^'(r), and the brackets 
(...) indicate averaging over random scattering potential 
between the puddles at a given shape of the puddles. The 
only, but crucial difference with the conventional case of 
s-n junctions (See, for example, l^, 14 1), is the bound- 
ary conditions for Eq. |4] at the normal-superconductor 
surface, which determine the sign of r]i. 

For the case when the size of the puddle is larger 
than the coherence length and the Andreev reflection 
on the puddles is effective the boundary conditions for 
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on the d-n boundary have been derived in Ref. 
Since the relevant energy for computing the Joseph- 
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son couphng, e « Dtrl\Yi — r^l^, is much smaller than 
the value of the order parameter in the puddles, the 
boundary condition for 9{v, e) is independent of e and 
depends only on the angle between the unit vector par- 
allel to the direction of a gap node oa and a unit vector, 
n(r), normal to the boundary at point r at the surface, 
: 9s{e,v) = /[a(r)], sin[a(r)] = ri(r) • ha- Here f{a) 
is a smooth, approximately odd and periodic function, 
/(a) « —f{—a), f{a) ~ f{a + tt), which grows from 
/(a) « at a = 0, to f{a) « ±C for a — 7r/4, where 
'--^ 1. Solving Eq. [4] with these boundary conditions, 
and using the standard procedure of calculation of the 
Joshepson energy we get 
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and C ~ Geff^^: V is the puddle volume, the integral 
is taken over the surface of the ith puddle, and Geff is 
the conductance of a metal of a size of order of the size of 
the superconducting puddle. In this case the magnitude 
of the 5-component of the order parameter generated at 
the superconductor-normal metal boundary is of order of 
the magnitude of the d-wave component. Thus it is not 
surprising that the value of J^^ in Eq. [1] turns out to be 
of the same order as in the case of SNS junction. 

If the distribution function of the mean-free paths is 
unbounded, and with certain probabilities one can find 
arbitrary large values of l(r), the mean field supercon- 
ducting solution always exists. However, if the puddle 
concentration is small enough, the transition from the 
state with global s-wave symmetry to the normal metal 
is triggered by a competition between the inter-puddle 
Joshepson coupling energy and the thermal (or quantum) 
fluctuations. Thermal fluctuations destroy the coherence 
between two puddles when Jij kT, which gives us an 
expression for the critical temperature Tc of the s-wave 
superconductor-metal transition 



T 91. 

where R is the inter-puddle distance. 
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We would like to stress that the existence of the s-wave 
superconducting phase is a generic property of the sys- 
tem because the long-range nature of the decay of Eq. [5] 
ensures that near the superconductor-normal metal tran- 
sition and at small enough temperatures the supercon- 
ducting puddles are separated by a distance larger than 
their size. 




S-wave\ normal metal 
disorder 



FIG. 1: Schematic phase diagram for the case when d-wave 
superconducting state is destroyed as a function of increasing 
disorder strength. 



In principle, the situation described above can be real- 
ized when grains of rf-wave superconductors are embed- 
ded into a normal metal artificially. In random systems 
the critical point can be identified by finding the set of 
"optimal puddles" which lie on the critical links of "the 
percolating cluster". In this case the properties of the 
s-wave phase and the dependence of the critical temper- 
ature Tc on the parameters of the system depends on de- 
tails of the distribution function of the disordered poten- 
tial. To illustrate the situation we consider here a simple 
model where the mean free part Z(r) is a random function 
of coordinates with a Gaussian distribution characterized 
by an average a variance ctZq, and a correlation length 
which is of order ^o- To be concrete, we consider the 2D 
case. Then the distance between the puddles becomes of 
order of their size, the amplitude of fluctuations of the 
order parameter becomes of order of the average, and the 
system has a transition to the s-wave state when I ~ Id 
and T < Td 
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Here Tco is the critical temperature of a pure d-wave su- 
perconductor, li Iq — I ^ la^ the distance between "the 
optimal puddles" is much bigger than their size. We can 
characterize such puddles by a value of the mean free 
path lopt > lo averaged over the volume of the pud- 
dle. In this case Aopt ~ Aolo/{lopt - hY^"^ < Aq, 
the size of the puddle is of order of the zero temper- 
ature coherence length ^opt ~ £,^1 1 {lopt - loY^'^ > ^o, 
and the characteristic distance between the puddles is 
of order of Copt exp[(?opt - /2a'^T{lopt - k)]- Here Aq 
is the magnitude of the order parameter in a pure d- 
wave superconductor at T = 0. This expression has 
a minimum at {lopt — lo) {lo ~ Oi ^^'^ therefore 
Ropt ~ exp[{lo - l)/loa'^]. Using Eq.glwe get 



Tc ^ rcoCTexp[- 



{lo-l) . 
lo'J^ ^ 
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At very small values of Tc the phase transition between 
the s-wave superconducting phase and the normal metal 
is triggered by quantum fluctuations of the order param- 
eter. In this case to determine the point of the transition 
one has to compare Jij with the " quantum temperature" 
characterized by the zero temperature superconducting 
susceptibility Xi of individual puddles [l^, [itI . [isj . Their 
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ing. There are a number of interesting predictions we 
can make. 1) There should be a transition from a glob- 
ally d-wave to a globally s-wave superconducting state at 
a doping concentration x — X2. (Some evidence of such a 
transition may already be present in the experiments of 
Ref. 1^.) 2) In the metallic state with x > X2, the conduc- 
tivity at low temperature should diverge as x X2, the 
Hall resistance should vanish, and the Weideman-Franz 
law should be increasingly strongly violated. 
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D. Fisher, A. Kapitulnik, and P. Young. The research 
was supported by NSF Grant No. DMR-0704151, and 
DOE Grant DE-AC02-76SF00515. 



FIG. 2: A qualitative illustration of the global d-wave to s- 
wave transition. Solid lines represent boundaries of d-wave su- 
perconducting puddles embedded into a normal metal. Hatch 
marked areas indicate the areas were the s-wave component 
of the anomalous Green function -F°(r, r) is positive. Outside 
these areas F''{r,r) is negative, a) The case of small puddle 
concentration when the system has s-wave global symmetry, 
b) The case of big puddle concentration when the system has 
a global d-wave symmetry. 



finite value is associated with exponentially rare events 
of tunneling of the superconducting order. Consequently, 
the value of Xi turns out to be exponentially big while 
the Joshepson coupling between the puddles decays with 
the inter-puddle distance only as a power law. There- 
fore, a generic feature of the s-n transition is that it takes 
place when the distance between optimal puddles is ex- 
ponentially big 18[. As a result Eq. [8] holds down to 



exponentially low temperatures. 

The cuprate high-temperature superconductors are the 
best established example of a d-wave superconductor. 
Here, the critical temperature, Tc, is known to vary 
strongly as a function of the doped hole concentration, 
X, producing two quantum critical points at which Tc 
vanishes: a lower critical doping concentration, xi, on 
the "underdoped" side, and an upper critical concentra- 
tion, X2 , on the "overdoped" side of the phase diagram. 
On the underdoped side of the superconducting dome, 
with increasing underdoping, these materials frequently 
appear to undergo a superconductor to insulator transi- 
tion @, 0, S]. Thus, the present considerations are not 
applicable. We assume, some of the more robust of our 
findings apply to the cuprates as Tc ^ with overdop- 
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